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Abstract 

We develop a finite temperature mean field theory in the path integral 
picture for an extremely dilute system of interacting Fermions in a plane. 
In the limit of short ranged interaction, the system is shown to undergo a 
phase transition to a superconducting regime. Unlike the well known BCS 
transition for metals, this phase transition is found to be very sensitive 
to the two dimensional nature of the problem. A quantitative estimate 
of this sensitivity is carried out by repeating the analysis at non - zero 
thickness. The validity of the mean field results are also proved using 
renormalization group techniques. 



1 Introduction 

It is well known from the theory of superconductivity^, that the low energy 
effective theory for electrons in metals is susceptible to a BCS phase transition 
at T = 0. The low energy theory which involves only the modes close to the 
Fermi surface (Kp) turns out (in any dimension D > 2), to be an effective 
one dimensional theory that describes a large number (N) species of Fermions 
satisfying a linear energy momentum dispersion relation. The inter - electron 
interaction for this theory turns out to be short ranged and attractive [|. The 
one dimensional nature of the effective theory and the linear energy momentum 
dispersion relation, are key ingredients that make the theory asymptotically 
free, while the attractive nature of the inter - electron potential is responsible for 
the formation of the gap in the energy spectrum. The effective one dimensional 
nature of the BCS theory is related to the fact that in considering modes close to 
the Fermi sea, (a natural cut off for the theory), one finds that all the directions 
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in momentum space tangential to Kp are degenerate, in the sense that the 
energy depends only on the directions normal to the Fermi sea. Hence one 
naturally obtains a large degeneracy, denoted by the number N, which turns out 
to be the area of the Fermi surface measured in units of the cut off. Moreover, 
this dimensional reduction to one dimension can be carried out in a manner 
which is insensitive to the true dimensionality of the problem which is one of 
the reasons behind the universality of the BCS transition]^, [I], 

The attractive inter - electron potential is also generic to all metallic systems 
in D > 2. This can be best understood by appealing to renormalization group 
(RG) arguments. RG analysis shows that due to coulomb screening all long 
ranged interactions (repulsive or otherwise) behave in a screened or short ranged 
manner. If one now decomposes the short ranged interaction into its various 
angular momentum components, then the RG flows for the various components 
couple to each other. The structure of these coupled equations suggests that it 
is always possible to find some value of I (the angular momentum ), for which 
the effective coupling for that sector will flow to negative values. The couplings 
that do run negative , i.e become attractive, lead to a BCS transition, while the 
others renormalize to zero. This is the so called Kohn - Luttinger effect j|, |(j|. f7|. 
The Kohn - Luttinger effect, and the dimensional reduction to one dimension 
make the BCS phenomenon truly universal, irrespective of the true physical 
dimensionality of the metallic systems or the microscopic details of the inter 
electron interaction. 

Recently there has been much theoretical work towards understanding the 
behavior of two dimensional dilute electron - hole systems^], [|. Whether or 
not such systems can undergo a phase transition to a conducting or supercon- 
ducting regime, and the dependence on the true physical dimensions of such 
possible phase transitions is the object of much speculation, as these systems 
arise naturally in the context of high T c superconductors. The physical simplifi- 
cations mentioned before in the context of the usual BCS theory do not apply to 
these dilute systems, where the physical momentum scale A >> K F . For such 
systems it does not make sense to linearize the energy - momentum dispersion 
relation around Kp and reduce to a one dimensional problem. We present here 
the analytical study of such a dilute system of electrons, interacting through an 
attractive short ranged potential. We carry out the study at finite temperature 
in the absence of impurities or disorder, and find that such systems undergo a 
phase transition to a superconducting state for arbitrarily small values of the 
attractive coupling. Moreover we find that they do so in a manner which is 
extremely sensitive to the two dimensionality of the problem. 

We will start with a short review of the quantum mechanical two body 
problem involving an attractive delta function potential in two dimensions, as 
it is an excellent toy model that encapsulates all the essential features of the 
corresponding many body problem, eg. the formation of bound states and 
asymptotic freedom. We shall then focus on the second quantized version of the 
problem, which will be analyzed in the path integral formalism. This will be 
followed by an exploration of the sensitivity of the BCS like phase transition to 
two dimensionality. In a separate section we shall justify the use of mean field 
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theory in the context of this problem using renormalization group techniques. 

2 Short Ranged Attraction in two Space Dimen- 
sions and Asymptotic Freedom 

The Two body problem: 

Let us consider a system of two particles interacting through an attractive short 
ranged potential in two dimensions ju| [ll], This two body system has an 
ultraviolet divergence, and the subtleties associated with the renormalization 
of theories which have one marginal parameter. The Schroedinger equation for 
this system (in momentum space) is: 

(f^ - #)*A(P) - ff(A)M(p) / PA(q)^A(q)[d 2 q] = (1) 
Note: [d 2 q] is a short hand for |^|. 

The equation above is regularized by the presence of pa(p) = ©(M < A). In 
the absence of the regularisation, the equation is ill defined, because the ground 
state energy of the system is not bounded below. This is an artifact of the 
scale invariance of the equation which implies that if E is an eigenvalue, then 
so is sE, where s is a scaling parameter. Hence the bound state energy can be 
either 0, or — oo. We can make sense of this scenario, by formulating the basic 
dynamical equation with a cut-off, as above, and ask the question as to how 
the dimensionless coupling constant g depend on A, such that all the physical 
quantities have a finite limit as A — > oo. Solving for ^a(p) gives: 

*a(p) = -^-=pa(p) (2) 

n E 

Where A a = g{A) J pa(p)^ A{p)[d 2 p]- Putting the expression for \&a(p) m 
the equation for A a, we have: 

g-HV-f #%<^] = (3) 

2m 

This expression for g~ 1 (A) is log-divergent. The situation can be remedied if 
we choose g(A) in the following way: 

<r x (A)= / pioo -g- 1 ^ [d 2 P ] (4) 

J 2m + 2m 

With this choice of g, Eqn(3) above has a finite limit as A — > oo, and the 
resulting theory is asymptotically free. The essence of renormalization here is 
the trading of the bare coupling constant g in favor of the number p, which 
has the dimensions of momentum. This number is the true physical parameter 
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of the theory, which describes the strength of the interaction between the two 
particles. 

The nature of the divergence, does not change if we make the transition 
from the two body problem, to the corresponding many body problem of non- 
relativistic particles in a plane. Indeed, it can be shown, that the divergence 
present at the level of the mean field theory of the many body problem is the 
only divergence that requires renormalization. In other words, renormalizing 
the mean field theory is enough to cast the problem in a finite form. Moreover 
renormalization group arguments (which we present in the last section)show 
that the predictions of the mean filed analysis for this problem are qualitatively 
correct. In the next section we will generalize the two body Hamiltonian given 
above to the corresponding many body case, and recast the problem of evaluat- 
ing the thermodynamic partition function for this Hamiltonian in the language 
of path integrals. In the next section we shall use this path integral point of 
view to construct the mean field theory for the problem. 

The many body problem: 
The generalization of the Hamiltonian given in Equation (1) to a second quan- 
tized language is, 

H = J [rf D r]*^( r )[-V 2 ]* ai (r) - a(A) j [d D r]^ (r)*+° a (r)*„ 2 (r)* ai (r) 

(5) 

We are keeping the dimensionality of space (D) arbitrary for the moment. This 
Hamiltonian can be written in a manner more suited to a mean field analysis, 
by using the following two component spinors, 



(7) 



In terms of the Hamiltonian can be written as 



H 



= /d D (r)[* t (-V a ) 0S *-a*tf ^ )* 
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The composite operator ipt(r)ip^(r), and it's Hermitian conjugate appearing in 
the Hamiltonian above are the creation and annihilation operators for the quasi 
particles of the theory, which are scalars in the present case. The interesting 
physical features of the theory, such as the formation of a gap in the spectrum 
are related to the fact that (V'jV'j) 7^ ° m the ground state. Hence we would 
like to reformulate the problem in a language where these composite operators 
are the chief dynamical variables. We shall do it below in the finite temperature 
limit and in the path integral picture. 

The finite temperature features of the theory can be probed by calculating 
the thermodynamic partition function Tr[e~@ H ] for the Hamiltonian H given 
above in (5). Now for any normal ordered Fermionic Hamiltonian one can 
reduce the problem of calculating the partition function to the evaluation of a 
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Path integral over Grassmann variables, using the prescription, 

Z = Tr\e- !m \ = J e- Io dt J d D r[^(rt)^(rt)-H(^(rt)Mrt))] £[^t]£[^] (g) 

Here it is understood, that the arguments of H appearing in the path integral are 
the Grassmann variables corresponding to the Fermion creation and annihilation 
operators appearing in Equation (5). This functional integral can be rewritten 
in the following way. 

Z= e ° V P J D[¥]D[*]D[^}D[<P} 

(9) 

The functional appearing in the exponential is to be thought of as the action S 
of an interacting theory. Here </> is an auxiliary complex scalar field. To see that 
this is the correct action, we can integrate out the scalar fields, and recover an 
effective action for the Fermi fields which is the action appearing in (8). i.e. 

J D[cj)]D[^]e- s = e - SF '^ (10) 

where 

S Fermi = P fdt J d D r[^\r)(^d t -W 2 )a^(r)-2g 2 r\r)^(r)Mr)Mr) 

(11) 

It is now obvious that S Fermi is the action generated by the Hamiltonian in (5), 
when translated to the path integral picture using (8) if we identify 2g 2 with 
the coupling constant a(A) appearing in (5). A comparison between equations 
(7) and (9) makes it obvious that (j> describes the quasi bosons of the theory, 
i.e (j) = V'lV'T- Reformulation of the theory in terms of the composite operator 
describing the quasi scalars can now be accomplished, by integrating out the 
Fermi fields appearing in (9). This produces the effective action for the complex 
scalar field. 



S Efs [4>,4> ] ,P]=p\ J d D r\0\ 2 -Trln 
Transforming to Fourier space, we have; 



(^-V> 3 







(12) 



(l 9t _V> 3 W f )' J' (2^+1)^2 +/3^ 

(13) 

The product on the left hand side is over the (odd) Matsubara frequencies. The 
evaluation of this trace can be carried out as follows; let 

A = lnn»(l+ .^a, 2 ) (14) 
(2n + \y-K z + y z 
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Where x = f3g\<j>\ = f3A, and y = (3k 2 . Now; 



- A — y°° i _ nycc \ /-i r\ 

dx 2 -°°(2n+l)% 2 + [3 2 E 2 1 (2n-l) 2 TT 2 +[3 2 E 2 { ' 

In the equation above, E = (k 4 + A 2 ) 1 / 2 . Now recalling the identity 

lanh( ») . if Sri _J F _ J (!„ 

we get 

^A=^tanh(^) (17) 

, 2 I0I 2 tanh(M) 

V 2 ' dx 2 (18) 
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E(x) 

Where S(x) = + w 2 (fc)) 1/2 , u(k) = k 2 . 

Mean field theory: Mean field theory corresponds to the saddle point of 
the non - linear scalar field theory described by (12). To get the saddle point 
equations, it is sufficient to consider a constant value for the scalar field <j> c . 
Using the form given in (18) for the functional determinant, the effective action 
for a constant value of the scalar field is, 

Where F denotes the free energy of the system. 
The mean field (saddle point equations): 

The equation for the saddle point is the equation for the extremum of the 
effective action, or the Free energy, hence we have from the form of Se// above, 

SF _ , 2 f , D tanh( ^ 2 +^» 1/2 ) 



-Q^!-; 2 [ d °! ^ y 2 ) (20) 



This is the familiar gap equation, and A = \g\<j> c has the physical interpretation 
of the gap in the energy spectrum. This equation is divergent in any dimension 
D > 2. The divergence can be made explicit by considering the zero temperature 
(/? — > oo ) limit, where it becomes, 



"7 <■*(»))■/» -' (21 » 

In the case of D = 2, which is the critical dimension for the thcorythis 
ultraviolet divergence is logarithmic, and it can be renormalized if we let g~ 2 ~ 

A 2 

f Q Ji , where is a dimensional parameter that sets the scale for the 
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rcnormalized theory. After incorporating this renormalization, equation (20) 
above in (2+1) dimensions cam be written in a manifestly finite form.; 




tanh(/Vu 2 + A 2 /2) - 1 
(a; 2 + A 2 ) 1/2 



= ln( — ) = Lim\2 




v/^ 2 + M 2 Jo V^ 2 + A 2 
(22) 




Here lu = k 2 . In terms of the dimensional variables y = — and £ = /3A, the 
renormalized gap equation above becomes; 



This is to be thought of as the equation of state for our system. The other 
dimensionless quantities like the temperature in the units of the binding energy 
x = flu = ^ can be recovered from the parametric equation above. This equa- 
tion is the same as the familiar gap equation for the BCS superconductor, and 
the numerical values for the physical quantities e.g. kT c /A predicted by the 
above equation of state are exactly equal to those of the usual BCS supercon- 
ductor. Hence we observe, that in the absence of impurities, a dilute system of 
electrons, interacting through arbitrarily small, short - ranged attractive poten- 
tials, will undergo a phase transition to a superconducting stage. 

Although the equations describing the phase transition are very reminiscent 
of the BCS transition for metals, they do in fact describe a different physical 
situation. This is evident if we note that the arguments that led to the mean 
field equation for this dilute system of electrons are very special to the two 
dimensional nature of the problem. Indeed the nature of the divergence of the 
gap equation itself is different, in higher dimensions, where the equation does 
not admit any solutions for values of the coupling constant less which are less 
than some critical value. 

This is different from what happens in the theory of metals, at the BCS 
transition point. As mentioned before, in the case of metals, one deals with 
an effective theory of modes near the Fermi surface, which can be reduced to 
a one dimensional theory with a linear dispersion relation which leads to a 
similar logarithmic divergence. Moreover, this reduction is insensitive to the 
true dimensionality of the problem (the exceptional case being D = 1) which 
makes the physics of the phase transition blind to the true dimensionality of the 
metal being studied. 

It is clear from the discussion so far, that the model we are studying now, 
does not admit of these simplifications, because of the diluteness of the system. 
Since for our system, ~ 0, where Kp is the Fermi momentum, and A is the 
typical energy scale of the problem, it is not meaningful to look for a theory 
for the modes close to the Fermi surface. This in turn ( as we pointed out 
above ) makes the phase transition depend critically on the two dimensionality 
of the model, which is an idealization. So we now probe the effects of a finite 
transverse spatial direction on the model. 
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3 The Analysis at Finite Thickness 



In this section we ask the question as to how thin the system has to be for it to 
be considered two dimensional. To keep matters simple, let us investigate the 
situation at zero temperature. In particular we want to see the transition from 
the three dimensional case to the two dimensional one. Eqn (21) tells us that in 
dimensions greater than two, the strong and weak coupling phases of the system 
do not match smoothly, i.e. there is a critical coupling g~ 2 r ~ A D ~ 2 , below 
which the gap equation does not admit any solutions. This is indicative of a 
phase separation between the strong and weak coupling regimes. To understand 
the effects of a finite thickness L3 on the system, we will address the following 
questions. 

a: What is the correct renormalization prescription to use for g~2, for small 
but non zero values of the the thickness L3. 

b: Given a finite thickness of the system L3, and the fact that the system is in 
the strongly coupled phase, how should L 3 approach zero, so that the gap (A) 
remains finite. 

At finite thickness, the non-relativistic dispersion relation reads u>(k) — k 2 + 
(tt) 2 ' wnere n takes on integer values. So the gap equation at finite thickness 
is; 

^I^^wTiEff^Twr 1 (24) 

fco's are the Matsubara frequencies which take on continuous values in (—00, 00) 
at T = 0. 

Summation of the series: 

The series given by S = J2„ -/(«*>„), where f(iun) = (fc 2_ (l ^)2 +fc 2 +A 2 
and oj n = (^p), can be evaluated by considering the following contour integral; 

If the contour of integration is taken to be the infinite circle centered at the 
origin, then the integral vanishes, i.e. Res(f (z)n(z)) = 0. The poles of the 
integrand are at z = ioj n , and at the four roots of ((fc 2 — z 2 ) 2 + k 2 + A 2 )). Let 
us denote the residues at the last four poles by Ri. The residues at the poles 
labelled by the integer n are j^f(iu n ). Hence S — L^^ 1 Ri. Summing the 
series by computing the four residues, the gap equation at finite thickness reads 
as; 



coth[(a - ib)^ 2 L 3 /2} coth[(a + i6) 1 / 2 L 3 /2] 



(a-ib) 1 / 2 (a + ib) 1 / 2 



(26) 



Here a = k 2 and b = (fc 2 , + A 2 ) 1 ^. 

The question now is to chose a suitable renormalization prescription for the 
coupling constant. To do this we must first understand the nature of the diver- 
gence in the equation. As coth(a;) ~ 1, for x » 1, equation (26) tells us that 



8 



the nature of the divergence is linear for any finite value of L 3 and logarithmic 
for L 3 = 0. Hence at a finite thickness, choosing g~ 2 ~ J A [d 2 k] ( fc2 + *2)i/2 , i- e 
trading the dimensionlcss coupling constant for the dimensional parameter /x 
should make the problem manifestly finite. 

To carry out the renormalization, we shall isolate the divergent part of the 
integral appearing in the gap equation (27). Recalling that the divergence in 
the integral remains the same when we let L 3 — > 00, i.e. let the hyperbolic 
cotangent go to 1, we can add and subtract, this divergent amount from the 
integral and rewrite it as; 
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[d 2 k][dk ]L 3 



CO 



th[(a-ib) 1 / 2 L 3 /2} 



+ 



1 



Anib 
[d 2 k}[dk ]L 3 
iirib 



(a - ib) 1 / 2 
1 



{a-ib) 1 ' 2 



(a - ib) 1 / 2 
- ((a - ib) -> (a + ib)) 



- ((a-ib) -> (a + ib)) 



Carrying out the integration over momenta (which are now finite), we have; 



-(i6) 1/2 L 3 



-(-i6)!/2L 3 



+ L 3 



[dk ] 

b l/2 



(28) 



The first integral in the equation above is finite while the second one encodes 
the linear divergence of the problem. It diverges as \fk~o, and since fco has 
the dimensions of fc 2 , it's divergence is linear in the momentum cutoff. So the 
appropriate renormalization prescription to use in this problem would be to let 

g~ 2 = L 3 J Q A ( fc 2_^2)i/4 ■ So the rcnormalized gap equation now is. 



4k/x/( A2 /M) - /( A2 /A)] = G(L' 3 , A 2 /A) 



(29) 



+ 
(27) 



Here f(a) = j Q a j^^, L> 3 - VAL 3 , and G(x, A 2 /A) = f< 
The above equation tells us that: 



A /A dk 



Arg 



-V»(feg+i) 1/2 » 

-^/-i(fc^ + l) 1/2 ^ 



G(^,A 2 /A) 
~TlJW7W 



/(A 2 / A) 

1WM 



(30) 



In the limit A — > 00, G remains finite, and the ratio of the two linearly divergent 
integrals approaches unity, and we are left with a result that is insensitive to the 
finite size, i.e. A = fi. So we see that for any finite value of L 3 , the behavior of 
the system is the same as that of a truly three dimensional system in the strongly 
coupled phase. This is indeed signalled by the renormalization prescription used 
above i.e g~ 2 <~ A, which is precisely of the order of the critical coupling in the 
full three dimensional theory. 
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The interesting dependence on the finite size is in the opposite direction, i.e. 
in the approach towards two dimensionality. So let us now ask the question, 

how should L 3 approach zero (when the cutoff is removed) so that remains 

finite ? To answer this question, we revert back to Equation (26), and change 
the renormalization prescription for g. We let g be renormalized in the way it 
would be for a truly two dimensional theory, i.e. g~ 2 = ln(^j-), and estimate 
the allowed range of values for L 3 , which is now thought of as a function of the 
cut-off. Expanding the hyperbolic cotangents in equation (26) to first order in 
L 3 , we get, 

1 OA2 r 11/" 1 

_!„,_, . jy k]mwTW _ + _ Ll jy A]m __ (31 ) 

or, in terms of dimensionless quantities, 

InA = CL\t? / -™= (32) 
M Jo V x 2 + 1 

Here C denotes a positive number (whose precise value is unimportant). Hence 
from (31) above, we have, 

hi(-) 

AL 3 = . V; (33) 
^ln(f) 

The equation above has the expected behavior, i.e, L3 — > 0, as A — > 00. A -1 can 
be though of as the measure of the lattice spacing in two dimensions. Hence 
we see that the theory will behave like a two dimensional theory even for a 
finite thickness, as long as the transverse thickness, L3 measured in units of the 
lattice spacing in the two dimensions is small. An estimate of this smallness is 
provided in the equation above. 



4 Justification of the mean field theory, a renor- 
malization group analysis 

In the previous sections we carried out the mean field analysis of a system of non 
- relativistic Fermions interacting through a short ranged attractive potential, 
by constructing the corresponding Landau - Ginsburg theory. At the level of the 
mean field theory, we identified a logarithmic divergence in the problem, which 
was removed by a renormalization of the coupling constant. The mean field 
analysis also predicted the existence of a gap in the spectrum, for arbitrarily 
small values of the attractive coupling. So the natural question that arises is that 
how can one justify the validity of the mean field approximation, which leads to 
these non - trivial consequences? In other words, we know of several examples 
where the predictions of mean field theory are qualitatively wrong; the D = 1, 
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Heisenberg chain at half filling^ is a good example of this; examples such as 
these provide us with a warning that motivate a surer justification of the mean 
field analysis. In this section here, we shall carry out a one loop renormalization 
group analysis of the theory, and show that the BCS instability predicted by the 
mean field analysis is also implied by the renormalization group flow equations 
, and this lends substance to the qualitative nature of the predictions of the 
mean field theory. 

The logistics for carrying out the renormalization group analysis for our 
model are going to be as follows; we shall start with a two dimensional, free, 
scale invariant theory, which we will define with a fundamental cut-off A. The 
Gaussian action is; 

So = J [dk ] J [d9] k[dk][^(k k)(ik - k 2 )^(k k)} (34) 

The action is invariant under the scaling transformation k — > sk,ko — > 
s 2 k , * — > s -3 *. 

Next we shall introduce various scale invariant perturbations around the free 
action, and study the evolutions of these perturbations under the RG transfor- 
mations. To generate the RG transformations, we shall split up the phase space 
(the k space ) into 'high' and ' low' degrees of freedom, the high and low degrees 
of freedom being defined as, ^high — [^(k), for (A/s < k < A) and otherwise], 
^low = \}&{k), for (0 < k < A/s) and otherwise ]. s is a real positive number 
greater than one. Following this we'll proceed to integrate out the high degrees 
of freedom, and get an effective theory for the low modes. This process will 
generate the RG flow equations for the various couplings that perturb the free 
action. Of special interest to us is the flow corresponding to the four point 
interaction, as that encodes the possibility or (lack of) a BCS type of instability 
in the system. 

Before we proceed to carry out the renormalization group transformation 
for the four point function, it is worth pointing out that unlike the usual BCS 
theory, we are not building an effective theory for the modes close to the Fermi 
surface, indeed we are interested in the scenario, where the system is sufficiently 
dilute such that ^f- « 1. In the case of the BCS theory, since the physically 
interesting degrees of freedom are the modes close to the Fermi surface, one 
is justified in linearizing the energy momentum dispersion relation around the 
Fermi energy and keeping the contribution only from the direction normal to 
the Fermi surface, as these correspond the degrees of freedom of significance 
to the RG transformations. The degeneracy of the angular directions, and the 
linearization of the dispersion relation allows one to reduce the theory to a 1 + 1 
dimensional theory with a relativistic dispersion relation, involving N species 
of Fermions, where N is a large number proportional to the area of the Fermi 
surface measured in units of the cutoff. These generic arguments carry through 
in any dimension greater than or equal to two, which is an explanation for the 
fact that all metallic systems, irrespective of the number of space dimensions 
face the BCS instability in the limit of T — > 0. 
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Figure 1: Graphs Contributing to one loop (3 function 



In the case we are interested in such simplifications are absent, although 
the BCS instability is present, however in a way, which is special to the two 
dimensional nature of the problem. 

To see this, let us start with a generic four point coupling; 

S int = Jn i [dk^[dk 0l }U(h)^\A)^H3)^(2)^(l)S(k)5(k ) (35) 

In the above interaction the number 1...4 denote the momenta and Matsub- 
ara frequencies. U is the four point coupling function, which in general will have 
a dependence on the momenta. The delta functions ensure than the momenta 
and the Matsubara frequencies corresponding to 1 and 2 equal those of 3 and 
4. All the momentum integrals run between and A. If we expand the four 
point function in a Taylor series, only the constant part (Uo) turns out to be 
of importance, as the terms with dependence on the momenta, correspond to 
irrelevant couplings, i.e, only the constant part leads to a four point interaction 
invariant under the scaling defined above. 

If we now split up the fields into the high and low modes, then the one loop 
corrections to the coupling constant are given by the three diagrams above (fig 
1), where the internal loop momenta correspond to the high degrees of freedom, 
i.e A/s < loop momenta < A. i.e. 



A 

A/s 



5Uo = ul1 dho1 } d2k] I d9[ {iko - (fc + Q?) + {iko - (k + cry) 



1 



2 (iko -k 2 )(-ik - {P-k) 2 y 

It is now obvious that the first two diagrams are identically equal to zero 
as the fco integrals vanish for these diagrams. Both the poles are on the same 
half of the complex ko plane because (k + Q) 2 and (k + Q') 2 being quadratic 
functions are always positive for all values of momentum transfer. Hence we can 
always close the contour in the other half plane and this makes these integrals 
vanish. The third integral however is non zero, and produces a flow, which is 
given by; 

/J=^ = ^ (37) 
dt 4tt v ; 



(36) 
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Here t = 



— Ins, and the solution to the flow equation is, 



Uo{t) = U °(°) (38 ) 

Thus we find that for an attractive microscopic coupling, (corresponding to 
Uo(0) < 0), a BCS instability is inevitable. Hence the divergence of the coupling 
constant that showed up in the mean field analysis is not an artifact of the mean 
field approximation, but is indeed a true signal of the onset of a BCS transition. 

We would like to emphasize the difference of this coupling constant renor- 
malization with its counterpart in the theory of metals. In metallic systems, 
the constraint that all momenta lie very close to the Fermi surface imposes a 
stringent dependence of the coupling constant on the transverse momenta, i.e on 
the directions orthogonal to the radial directions in momentum space. Moreover 
the only couplings that contribute to the flow are the ones for which k\ = —k 2 , 
and &3 = — &4, and Uo = Uo(kx.ks), i.e, it is a function of an angle. When 
this function is decomposed into it's various angular momentum components, 
the contribution from the first two diagrams (which are non zero, as the dis- 
persion relation is linear in the momentum) couple the flows corresponding to 
the various angular momentum sectors (the Kohn - Luttinger effect H, |(| Q]), 
in a manner that drive (for certain values of the angular momentum) couplings 
which might have been repulsive to negative values. These couplings then lead 
to the BCS instability according to the flow equation (37) given above. 

Hence it is clear that the system we are presently interested in is different 
from the usual BCS superconductor, although it also possesses a similar ground 
state. But there is no analogue of the Kohn Luttinger effect for this system, 
however a coupling that is attractive to start with will be driven by the RG flow 
to produce a BCS like ground state. Moreover, the arguments and the analysis 
we given above are very special to the two dimensional case as was shown in 
the previous section. 



5 Conclusion: 

We conclude that a dilute, two dimensional system of electrons, interacting 
through a short ranged attractive potential will undergo a phase transition to a 
superconducting regime. This phase transition is different from the usual BCS 
transition in several respects. 

a: The phenomena described here is very special to the two dimensionality 
of the problem. 

b: The phase transition is from a semiconducting / insulating phase to a 
superconducting regime. 

c: The diluteness of the system plays a very significant role in the analysis, 
as that is what warrants the use of the quadratic dispersion relation, for which 
the critical dimension is d = 2. 

On a speculative note, we would like to emphasize that these phenomena may 
be of relevance towards understanding the behavior of high T c superconductors 0, 
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In dimensions greater than two we find that this kind superconductivity will 
not arise for arbitrarily small values of the electron - electron coupling, and that 
the coupling will have to be of greater than some critical value g cr for the phase 
transition to take place. If one allows for a finite thickness of the material being 
studied, then we find that there are stringent bounds on the allowed range of 
thickness if the phenomena characteristic of two dimensions is to survive; in fact 
we find that the square of the thickness can exceed the lattice spacing in two 
dimensions only by logarithmic amounts. 
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